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SUMMART 


An approximate method of calculating pressures -within the tip 
region of a rectangular wing having a symmetrical circ ular -arc profile 
is developed hy means of -which it is possible to determine the pressures 
at supersonic speeds to a higher degree of accuracy than with the usual 
linearized methods. The method consists essentially of applying tip 
corrections realized from linear theory to the exact pressures calcu- 
lated for the two -dimensional flow region and then talcing into account 
the actual region Influenced "by the wing tip. A comparison of same 
calculated pressure distributions with experimental results Indicates 
good agreement for regions not affected by shock-boundary-layer 
interactions . 


INTRODUCTION 


Comparisons between theory and experiment indicate that, for 
rectangular wings of moderate thickness or at moderately high angles 
of attack, the pressure distributions calculated by linear theory are 
not in good agreement with experimental results. If viscous effects 
are neglected, the pressure distributions for the parts of the wings 
-where the flow is two-d imens ional can be calculated with near exactness 
by the method of characteristics. Other methods that provide less 
accuracy than the method of characteristics but greater accuracy +--h«u 
the linear theory are the shock-expansion technique (neglect of the 
reflections from initial shock that are considered in the character- 
istics method) and Busemann's second-order-approximation theory. For 
moderate Mach numbers and thiclmess ratios the shock-expansion method 
gives results which are nearly identical with those of the character- 
istics method and which are in very good agreement with experiment -vixen 
no strong viscous effects are present. For the region influenced by the 
wing tip, no theories of corresponding precision are as yet available. 
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In reference 1, Bonney introduces a method, of calculating the 
aerodynamic characteristics of rectangular •wings with an accuracy 
greater than. that possible with usual linear methods* The method, con- 
sists essentially of computing the pressures in the two-dimensional flow 
regions "by maans of Busemann *s second-order approximation and assuming 
that, in the region influenced "by the wing tip, the differences in 0 
pressure "between the wing surfaces decrease from those in the two- 
dimensional region to zero at the wing tip in accordance with the rate 
of loss in lif t found "by lin ear theory for an infinitely thin flat plate. 
The method is not suitable for the prediction of actual pressures in the 
wing tip region, however, because the increments in pressure due to wing 
thickness largely cancel and need not be computed far the tip region in 
order to determine the lift and moment characteristics of the wing. 

The present paper introduces an approximate , method by means of 
which it is possible to calculate the pressures in the tip region of a 
rectangular wing with a symmetrical circular -arc section to a higher 
degree of. accuracy than is possible with the usual linear methods. 

The method is basically similar to Bonney *s method, but the two- 
dimensional pressure distributions are calculated by more exact methods 
than Busemann. ’s second-order approximation, and tip corrections realized * 
from linear theory are applied to the pressure increments due to both 
angle of attack and airfoil thickness. In addition the actual flow 
field influenced by the wing tip is taken into account rather than the 
field between the free -stream Mach line and the wing tip as is done in 
linear theory and Bonney *s method. A brief comparison of same calcu- 
lated pressure distributions with experimental results for a 9-percent- 
thick wing at a Mach number of 1.62 is included. 


STMBOIS 

P pressure coefficient 

AP a pressure -coefficient 

AP 0 pressure -coefficient 

and shape 

p local static pressure 

p 0 free -stream static pressure 

cLq free -stream - dynamic pressure 

Xqp center of pressure, in fraction of chord from leading §dge 



increment due to angle of attack 
increment due to airfoil thickness 
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t/c 

c mc/2 

°n 

“c/2 


n 

M 

Mo 

a 


chord of airfoil 

maximum thickness ratio of airfoil 
section pitching-moment coefficient 

section normal-force coefficient 

section pitching moment about 50 -per cent -chord point, 
positive -when it tends to rotate leading edge of 
section upward 

section normal force, positive upward 

local Mach number 

free -stream Mach number 

angle of attack 




p = K 2 - 1 


x, y, z coordinates of mutually perpendicular system of axes 

£, tj coordinates •which replace x and y, respectively, used 

to indicate origin of source line 


u 


x-component of disturbance velocity, positive in flight 
direction 


dz/dx slope of airfoil surface 

Kj>, Kg, constants used to determine location of point in 
K4, K5 distorted flow field (defined in fig. 5) 


V 


free -stream velocity in flight direction 


Subscripts : 

1 pressure coefficient in region influenced by .wing tip 

determined by linear theory 

pressure coefficient in two-dimensional flow region 
determined by linear theory 
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k 


3 pressure coefficient In two-dimensional flow region 

determined by use of oblique -shock theory and Prandtl- 
Meyer’s equations for espansicn of two-dimensional 
supersonic flow 

n 

Subscript notation for u and — indicates the origin of source 

dx 

line in terms of x . and y, respectively, except -when noted otherwise . 

Subscript b on x and y in equation (5) indicates location of 
boundary between two-dimensional and three-dimensional flow regions. 

Primed symbols : 

Primes indicate the location or coordinates of a point in the 
distorted flow field. 


ANALYSIS 


In the linearization of the three-dimensional equations of 
potential flow the following assumptions are made: the total-pressure 

losses across shocks are negligible, disturbances propagate on the air- 
foil along straight lines inclined at the Mach angle of the free stream, 
and no interaction exists between the thickness and angle -of -attack 
effects. In order to improve the accuracy of the pressure -distribution 
calculations for the tip region these simplifying assumptions must be 
modified to include conditions that are more nearly in, agreement with 
actual flow conditions. 

The effects of the total-pressure losses across the leading-edge 
shock and of the interaction between the pressures due to thickness and 
angle of attack are approximated by the use of pressures that are calcu- 
lated for the two -dimensional -flow region by methods involving oblique- 
shock theory and the Prandtl-Meyer equations for the expansion of a two- 
dimensional supersonic flow. If the method of superposition of linear 
theory is followed, the two-dimensional pressures are divided into two 
increments. (See fig. 1.) Thus, 

P3 = £Pct 3 + APog (1) 

where the pressure increment due to thickness AP ff is taken as constant 
for a given airfoil and the effects due to interaction are concentrated 
in the pressure increment due to angle of attack . AP a . The pressure 
increment due to interaction between angle -of -attack and thickness 
effects may be defined as the difference in pressure found when the 
pressure at any point on the airfoil is computed by two methods. In the 
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first method, T&ich gives results that are in very good agreement with 
experiment, the pressure is computed "by the usual shock-expansion 
te chn ique in which the thickness and angle -of -attack effects are con- 
sidered together. In the second method, the pressure is computed "by the 
shock-expansion theory, hut the pressure increment due to thickness is 
calculated separately at a = 0° and to this increment is added the 
increment in pressure calculated for a thin flat plate at the appropriate 
wngiA of attack. The assumption that the interaction effect is concen- 
trated in the increment in pressure due to angle of attack was made in 
order to reduce slightly the amount of required calculations. Actually, 
some calculations ' indicated that a more reasonable assumption, such as 
one in which the interaction effect was equally divided between the 
increments in pressure due to thickness and angle of attack, would 
generally lead to results in slightly better agreement with experiment. 
The act ual differences in pressure resulting from the use of the two 
methods of palculation are usually negligible as regards the calcula- 
tion of pressures in the tip region.. 

According to linear theory, the ratio of the pressure in the tip 
region of a rectangular wing with a symmetrical circular-arc profile 
to the pressure at the corresponding chordwise station in the two- 
dimensional flow region is given, for a = 0°, by the expression 

- ~(jt - COS " 1 pi) + - i (cosh -1 ~-r + ^ i(jt - COS -1 pi j] 

APtf-L *P\ 1/ * C [ • p|y| p y\ xj_\ 

^ =1+ K—d 


(2) 

where the equation is valid for -1 S pi < 1. (See appendix: for 

x 

derivation of equation and fig. 2 far definition of pertinent symbols.) 
The corresponding pressure ratio due solely to the angle of attack of a 
flat plate of rectangular plan form (reference 2) is 

~rr~ = - COS -1 (l + 2pi) (3) 

APog jt \ xj 


where 


-1 < pi <1 
~ x 

When the method of superposition of linear theory is followed but the 
two-dimensional pressure coefficient calculated by rigorous methods as 
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a 1)886 is used., the theoretical pressure coefficient anywhere in the tip 
region is defined "by 


P 


" ^3 


£®«1 

APa 2 ^ a 3 APag 


( 4 ) 


Equation (4) may he regarded as a correction of the pressure coefficient 
in the tip region Pi due to -wing thickness or angle of attack as cal- 
culated by linear theory by the corresponding ratio of the exact two- 

Po 

dimensional pressure to the Un sar '•hvo- fliTOft-ns ir ma-T pressure — It may 

P2 

be noted that equation (2) is indeterminate at ^ = 0 (at the •wing tip), 

c 

but the pressures at the station actually have a finite value -which 
varies along the chord. The equation also changes in sign from minus to 


plus infinity at the chordwise station = 0-50 (for the circ ular -arc 

profile), and thus a discontinuity results in the pressure distributions 
derived from equation (U) if the exact two-dimensional pressure coefficient 
due to thickness does not change in sign at exactly the Bame chordviBe 
point as in the linear theory. "When such pressure discontinuities do occur 
they are ascribed to inaccuracies of the method and are neglected in the 
‘fairing of the, pressure-distribution curves. Usually, the discontinui- 
ties are confined to a region less than 5 or 6 percent of the chord in 


length and located at the midchord point of the wing. For = -1 and 1, 
equation (2) reduces to 1 and 0, respectively. Equation (3), on the 


other hand, indicates a conical flow field with zero pressure ratio 
at the wing tip. 


AP. 


“1 


AP, 


°2 


Calculation of the distortion of the field influenced by the wing 
tip to allow for the propagation of disturbances along Mach lines 
dependent upon the local velocities over the wing is somewhat more 
difficult and can be accomplished only by approximation. Among the 
methods of distorting the tip field that were investigated are: a 

linear distortion in a spanwise direction, with no changes being assumed 
in pressure at the wing tip; a linear distortion in a chordwise direc- 
tion with the pressures at the trailing edge held constant; dis- 
tortions in both the spanwise and chordwise directions based an various 
concepts that disturbances propagate along curved Mach lines. Most of 
the techniques investigated gave results that were in fair agreement 
with experiment and thus indicated the purely arbitrary nature of the 
distortions, even though in same instances they appear to be linke d with 
the basic assumptions of some of the linearized theories. The best and 
most consistent agreement with experiment was obtained by transforming 
the field by the method described in the following paragraphs. 
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According to Jones' linear theory (see appendix), the pressure at 
any point L(x,y) (fig. 3) an an arbitrary spanvise station EF is 
determined hy the pressure due to the action. of the line sources and 
sinks representing the wing surface from the leading edge to the 
point L(x,y) and by the decrease in pressure due to the action of the 
line sources and sinks representing the tip and originating between A 
and G. The pressure disturbances are assumed to propagate along curved 
lines similar to AD and the pressure decrements resulting from the tip 
effects calculated from the sources and sinks from A to G are 
assumed to be felt at the reference station at the point L'(x*,y). By 
equation (4), the pressure at the point L(x,y) is corrected so that, at 
the same chordwise station, the pressure at the edge of the region 
influenced by the wing tip is equal to the exact two-dimensional 
pressure. The decrease in pressure due to the wing tip is taken as the 
difference between the exact two-dimensional pressure at the chordwise 
station corresponding to the point L(x,y) and the corrected pressure 
at the point L(x,y) . 

Inasmuch as the curved Mach line AD from the leading edge of the 
tip lies at the boundary of the two -dimensional flow region, itB loca- 
tion can be determined accurately from the exact two-dimensional calcu- 
lations from the following relationship derived by Erick and Boyd of 
the Ames Aeronautical Laboratory: 



where M 3 denotes the local two-dimensional Mach number at the 
point x/c. The curved Mach line GL' is assumed to be equal in slope 
and curvature to the corresponding chordwise section of AD. In other 
words, the curve GL' is found by translating AD toward the right 
until it passes through G. The determination of the coordinates of a 
large number of points L'(x',y) is simplified by first finding the 
locus of the points T(x,y) (fig. that corresponds to those 
points F '(x',y) which lie on the trailing edge of the wing DC. At 
each spanwise station EF the desired point is located graphically by 
tr anslat ing AD to the right until it passes through F or F ', thus 
locating G. From G a line is drawn parallel to the Mach line AB. 
The Intersection of this line GT with the station line EF gives the 
locus of the desired point. It should be noted that T(x,y) is the 
point required to determine the pressure at the trailing-edge 
point F'(x',y)j therefore, the procedure for locating the point is 
reversed from that for locating a transposed point L’(x',y). Typical 
loci of points T(x,y) are shown in figure 4. The ordinate of the 
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point L’(x%y) (fig. 5(a)) ia terms of the ordinate of L(x,y) and 
some graphically determined constants is 



(6) 


This procedure amounts to expanding the region HT linearly to N 'F ' . 
The pressure at the point L’Cx^y) is then found from ‘ 


■^at x’,y C^3)at x*,y (^3)at x,y + 


K) 


AP. 


at x,y AP 


fi 

°2 


^ a 3)at x. 


AP, 

y ^ 


ag 

°2 


(T) 


In order to determine the pressures at span-wise stations lying inboard 
of the point -where the line or loci of the points - T(x,y) intersects the 
■wing trailing edge (for example, fig. 4, lower Burface), a subterfuge 
must be used and pressure ratios must be computed for points behind the 
model. This computation is carried out by assuming that equations (2) 
and (3) are valid without change and by computing the exact two- 
dimensional pressure for the point as if the upper and lower surface 
extended without a break and with the same curvature beyond the actual 
wing trailing edge- 


■When a slightly greater discrepancy between the calculated and 
experimental pressures, probably not exceeding twice that resulting 
from the use of the proceeding method, can be tolerated (in general the 
changes in pressures and in integrated characteristics will be small), 
the calculations may be simplified by assuming that the points T(x,y) 
coincide with the trailing-edge points F ’(x',y) and no tralling-edge 
distortion is present. With equal, accuracy for the outboard stations 
and considerably better accuracy for stations approaching the two- 
dimensional region as compared to the previous simplified method, the 
calculations may be simplified by distorting the tip region linearly in 
a spanwlse direction with the tip pressures held constant. In this case 
the pressures are calculated directly from equation (4) in which the 
AP„_ AP. 


CT 1 


AP, 


ttT 

and — are computed for the point P(x,y) and are 


cr 2 


AP, 


og 


ratios 
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applied at the point P’(x,y ’) • The relation between the initial and 
transposed ordinates is 

xi„x5t i 

c 

•where the quantities and Ke 

Because the method of calculating pressures in -wing-tip regions 
that is presented in this paper requires the use of Mach lines, it is 
obviously restricted to cases in -which the flow over the airfoil is 
always supersonic. The technique can also be applied to swept untapered 
wings of constant section thickness ratio, but this application has not 
been carried out herein because, in most practical cases, disturbances 
from the apex of the swept leading edge will affeot part of the region 
influenced by. the wing tip. The pressures in the field influenced by 
the apex of the leading edge or on wings having other than circular-arc 
sections can similarly be determined. In all oases only the form of 
equation (2) and possibly equation (3) will change. It Bhould also be 
noted that the accuracy of the pressures approximated for the tip region 
depend directly upon the accuracy with -which the two-dimensional pres- 
sures can be determined. 


c K5 

are defined in figure 5(b) • 


COMPARISON 'WITH EXEERIMEHT 


A comparison between same pressure distributions calculated by the 
method of this paper with experimental pressure distributions is pre- 
sented in figure 6 for a rectangular wing having a symmetrical circular - 
arc profile, 9-percent thick. At the Mach number used for this compari- 
son (Mo = 1.62) the station ^ = -0.282 lies approximately in the 

c 

middle of the region influenced by the wing tip. The agreement between 
the calculated and experimental values is generally good for regions not 
affected by shock-boundary-layer interaction or by a Blight misallnement 
of the upper surface of the flap. At the higher angles of attack a slight 
discrepancy exists between the calculated and experimental pressures ‘ 
over the forward part of the wing theoretically not influenced by the ' 
wing tip. This discrepancy is believed to be due to the occurrence of 
wing twist during the tests that was not' accounted for by the experi- 
mental method of measuring angles of incidence • 

Additional calculated pressure distributions for other spanwise 
stations are presented in figures 7 and. 8. The variation with spanwise 
location and angle of attack of the aerodynamic coefficients obtained 
by mechanical integration of the pressure distribution curves of 
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figures 7 and 8 are shown in figure 9* For comparison the section 
coefficients determined by linear methods are included in the latter 

figure . At station Z = -0 .282, the values of <% and Cm calculated 

c 

by the present method for a = 4*55° (integrated from fig. 6) -were 0.174 
and 0 *028, respectively, and the calculated center of pressure was at 
the 34.0 -percent chord point. The corresponding experimental values 
of Op and Cm were 0.200 and 0.031 with the center of pressure at 
the 35 *5 -percent- chord point. A brief comparison of some of the aero- 
dynamic characteristics obtained by the present method with those 
determined by Bonney's technique for a wing of aspect ratio 2 indicated 
very good agreement in normal-force coefficients and in pitching-moment 
coefficient at low angles of attack. At higher angles of attack, the 
present calculations show a forward movement of the center of pressure 
with a that is not predicted by Bonney *s method. 


CONCDJDIN& REMARKS 


An approximate method has been derived for calculating pressures 
within the tip region of a rectangular wing having a circular-arc 
profile . Comparison of results obtained by this method and by linear 
theory with experimental results indicates that the approximate method 
is more accurate than the linear theory far predicting these pressures 
at supersonic speeds. 


Langley Aeronautical Laboratory 

national Advisory Committee for Aeronautics 

L angley Air Force Base, Ya., July 28, 1949 
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APPENDIX 

PRESSURE RADIOS DUE TO "WING- THICKNESS AND ANGIE OF ATTACK 


According to linear theory, the pressure coefficient at any point 
on a fixed, -wing is given "by 


P = AP ff + AP a (Al) 

where AP^ is the pressure on the surface, of the given wing at zero 
angle of attack and. AP a is the pressure on the surface of a flat plate 
of the given plan form at the angle of attack a. The equation can he 
rewritten for convenience in obtaining the derivations in the main body 
of this paper as 


P = AP, 


AP, 
ff 2 AP 


Si 

02 


+ AP, 


AP, 


o-l 


°2 AP, 


02 


(A2) 


where P denotes the pressure in the tip region as calculated by the 
method of this paper and the subscripts 1 and 2 indicate the pres- 
sures in the tip and -two-dimensional flow regions, respectively, as 


calculated by linear methods. 


The ratio 


AP, 


m 


a l 


02 


was first derived by 


Busemann (reference 2) by conical -flow methods and may be expressed as 


AP, 


AP, 


1 -l 

; = — OOS 

0£ « 


0- + <) 


(A3) 


where 


-1 ^ pZ. £ 1 

x 

For a wing of rectangular plan form and circular-arc section, the 
pressure due to thickness can be derived most conveniently by the use of 
semi -infinite line sources and sinks by the method of Jones (reference 3) • 
Figure 2 shows the location and positive directions of the system of 
axes used in the derivation. The semi-infinite line sources repre- 
senting the airfoil surface are assumed to originate at the spanwise 
station -h and extend to the right. The line sources of equal 
strength but opposite sign required to simulate the square wing tip are 
assumed to originate at station y = 0 and extend to the right. 

Station -h is taken sufficiently distant from station y = 0 so that 
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the Mach cone from the point x = 0,y = -h -will not overlap that from 
the -wing tip. From equation (A4) of reference 4 , the disturbance 
velocity on the -wing may be expressed by 

11 = ' u O, -h " ^0,-h - -1*0,0 + ^*0,0 (A4) 

•where the subscript notation indicates the origin of the source linn and 
the u-expressians are given by 

u. (x,y) = X. M fit - cos " 1 E i Z - Ujl 
M' itp dx[_ x - Fj 


- -? y - r (A ‘ 

+ i ^ jit - COS"- 1 - P - (y . ■ 1 l) . 

P y - til. x - IJ 

J 

■where |,tj represents the origin of the elementary source lines. For 
the circular-arc airfoil. 


= 2* 
c 


d g z _4 t 

cb^ " c 0 

Within the accuracy of the linear theory, the- pressure -coefficient 
ratio is 

-2ui 

APgi ~ ui 

APas " -2us = **2 (A6) 

V 

Substituting the appropriate parts of equation (Ah) into equation (A6) 
results in 


= 1 + 


-u 0,0 + ^“0,0 

%-h "^0,-h 
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■where the denominator of the fraction is restricted to the two- 
dimensional flow field outside and to the right of the Mach cane from, 
the point x = 0,y = -h and the numerator applies to the region within 
the Mach cone from the wing tip. Substituting equations (A5) into (A7) 
and simplifying give 






(AS) 


■where the equation is valid for 
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Figure 2.— Location arid positive direction of axes used In derivations of 

pressure ratios. 
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Figure 3 .— Sketch showing method of transposing a point in the flow field 
of linear theory to the distorted flow field which allows for local . 
f velocities on the airfoil. 
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Figure 4.—' Typical boundaries of region influenced by wing tip after 
distortion to allow for local velocities on airfoil. Symmetrical 

_i_ 

circular-arc airfoil; J = 0.09J M «» 1.62; a = 3*35* 
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(a) a m 0.55°. 


(h) a e 1.55°. 


Figure 6.— Comparison of calculated and experimental pressure distributions in the tip region of a 
symmetrical circular-arc airfoil. — = 0.09: — = - 0 . 282 : M = 1 . 62 . 
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Linear theory 



Distance from wing tip, y/c 

Figure 9.— Comparison of some aerodynamic characteristics of the tip 
region of a symmetrical circular-arc airfoil as determined, from 

linear theory and the approximate method. ^ = 0 . 09 . 
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